Dynamical tunneling in macroscopic systems 



I. Serban^'^ and F. K. Wilhelm^ 

'Department Physik, Arnold- Sommerfeld- Center for Theoretical Physics, and Center for Nanoscience, 
Ludwig-Maximilians-Universitat, Theresienstr. 37, 80333 Miinchen, Cermany 
^IQC and Dept. of Physics and Astronomy, University of Waterloo, 
200 University Ave W, Waterloo, ON, N2L 3G1, Canada 
(Dated: February 1, 2008) 

We investigate macroscopic dynamical quantum tunneling (MDQT) in the driven Duffing oscil- 
lator, charateristic for Josephson junction physics and nanomechanics. Under resonant conditions 
between stable coexisting states of such systems we calculate the tunneling rate. In macroscopic 
systems coupled to a heat bath, MDQT can be masked by driving-induced activation. We com- 
pare both processes, identify conditions under which tunneling can be detected with present day 
experimental means and suggest a protocol for its observation. 
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The phase space of a classical system can have forbid- 
den areas even in the absence of potential barriers, e.g. in 
the presence of external driving. Quantum-mechanically, 
these areas can be crossed in a process called dynamical 
tunneling [I] H]. So far, dynamical tunneling has been 
observed experimentally in microscopic systems, i.e. cold 
atoms [3] with very low dampig. Recent experimental 
progress has demonstrated many basic quantum features 
in macroscopic systems such as Josephson junctions or 
nanomechanical oscillators, overcoming the limitations 
posed by their coupling to the environment. Important 
for this success was the ability to reduce noise and cool 
to very low temperatures. 

In this paper we discuss the possibility of macroscopic 
dynamical tunneling (MDQT) i.e. involving a macro- 
scopic degree of freedom, like the phase difference across 
a driven Josephson junction. Classically, for certain pa- 
rameters, this system has two stable coexisting oscilla- 
tions with different amplitudes. This driven system will 
feel the influence of its dissipative environment strongly 
even at temperature T — Q. We demonstrate that un- 
der experimentally accessible conditions the tunneling 
between the two classical states can indeed occur and 
be singled out from the background of classical switch- 
ing events. We suggest an experiment where MDQT can 
be directly observed. Our result can be applied to verify 
quantum physics in systems with weak nonlinearity such 
as nanomechanical oscillators. Quantum tunneling it is 
also a potential dark count error process in the Joseph- 
son bifurcation amplifier. Here the classical switching 
between the two driving-induced, coexisting states in a 
Josephson junction was used for high resolution disper- 
sive qubit state detection |U|5Jinj|7]. 

Dynamical tunneling (in the absence of an environ- 
ment) has been studied using WKB in the parametric 
driven oscillator |S]. Activation rates in the presence of 
an environment have been studied in bistable systems 
[9l [To] . Dynamical tunneling with dissipation has been 
described numericaly [llj and multiphoton resonances 



have been studied perturbatively [12]. 

We study a harmonically driven Duffing oscillator, as 
an approximate description of a wide range of macro- 
scopic physical systems ranging Josephson junctions 
[H [T3] and nanomechanical oscillators [TU [TS]. The 
driven Duffing oscillator is described by the Hamiltonian 
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where F(t) — i^o('E"^* + (e~"^*) is the driving field with 
frequency v. For sub-resonant driving v < U, and below 
a critical driving strength Fq < Fc two classical oscil- 
latory states with different response amplitudes coexist. 
Considering a Josephson junction with capacitance C, 
critical current Ic and driving current amplitude / we 
can identify x as the phase difference across the junc- 
tion, m = {h/2e)^C, n = ^2eIc/{hC), Fq = hl/{2e) 
and 7 = mfl'^/ 24:. 

Following the Caldeira-Leggett approach, we assume 
an Ohmic environment and describe it as a bath of har- 
monic oscillators 
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with spectral density J{lo) ~ tt \f5{uj—uji)/{2miijji) = 
rnKw exp(— cj/wc) and lJc a high frequency cutoff. 

We transform this Hamiltonian using the the unitary 
operator U = exp{ii/t{a^ a + X^i ^I^i) similar to Ref. [S], 
where a and bi are the annihilation operators for the 
system and bath oscillators. Dropping the fast rotat- 
ing terms in the rotating wave approximation (RWA), 
we obtain 



2 



where, up to a constant we have 
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We thus obtain a time independent Hamiltonian at the 
expense of a form that is not separable in p and x. 
This transformation reduces the frequency f2 — fl6 
and increases the mass m — m/6 of the oscillators by 
6i = {uji — v) /uji in the case of the bath and S = 
+ Kcuc/ (Trfi^) for the main oscillator, where the 
additional term describes a deterministic force induced 
by dragging the system through its environment. 

We concentrate at first on quantum tunneling in the 
absence of bath fluctuations and study the system in the 
phase space. The classical Hamilton function Hq^\x,p) 
is portraited in Fig.[ljb) and (c) for a sub-critical driving 
strength Fq < = 2/9(2m^h^/"f)'^/'^. It has three ex- 
tremal points: saddle (s), minimum (m) and maximum 
(M) with coordinates {xc,Pc) in the phase space, where 
e € {m, s,M}. The curves satisfying Hq^\x,p) = E rep- 
resent classical trajectories. In the following we call E 
the quasi-energy. In the bistability region E E {E„-^,Es) 
where Eg = HQ^\xeTPo) there are always two periodic 
classical trajectories, around the two stable points (m) 
and (M) , with a small and large amplitude respectively. 

Using this phase space, we outline an experiment to 
observe MDQT during the transient evolution of the sys- 
tem. Without driving, the system relaxes to its ground 
state centered around (m). Then, after switching on the 
driving field one records the time needed for a transition 
to the large orbit as a function of a driving parameter 
e.g. frequency v. When two quantized levels pertain- 
ing to the two oscillatory states have almost the same 
quasi-energy, tunneling can occur, and enhance the total 
switching rate. 

We describe tunneling using the semiclassical WKB 
approximation which is an expansion in h close to the 
least action path. To find that path we solve the equation 
Hq^\x,p) = E and obtain four coexisting momentum 
branches ±pl,s{x, E) where 



Ps,Lix,E)=mh^ 



\ 37 



x^ T 



l8Fo 
37 



(4) 



with X = E/Fq — {7fLfl^)^/{6Foj). This configuration is 
reminiscent of Born-Oppenheimer surfaces in the molec- 
ular physics where dynamical tunneling has also been 
studied [T]. A real- valued psx corresponds to a classi- 
cally allowed area with an oscillating WKB wave func- 
tion, a complex-valued one to a classically forbidden area 
with a decaying wave function. At x ~ X, both trajecto- 
ries have the same momentum and position and connect. 
Here x = dpH'^\x,p) = but p ^ such that the 



the motion changes direction and continues on a differ- 
ent momentum branch. For all x < X both ps.l{x,E) 
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FIG. 1: Illustration of the Hamilton function and the "po- 
tential" landscape, (a) —p%^i^{x, E): the "potential" changes 
with E; classical turning points are found at p{xi,E) — 0. 
(b,c) h'^\x,p)\ in (b) white corresponds to high, black to 
low quasi-energy; (b): white lines corresponds to the L, black 
ones to the S branch; continuous lines correspond to real and 
dashed ones to imaginary valued momentum. 

are complex, thus this forbidden area does not infiu- 
ence the quantization rules within the WKB approxima- 
tion. The tunneling least-action trajectory which con- 
nects the two allowed regions only passes through the 
region x > X. Here the ps,L are either real or purely 
imaginary, i.e. —p^j^ S R. To study this region, we mir- 
ror the solution pl{x, E) around the X point as shown in 
Fig.[l|^a) and obtain a double well "potential" . The small 
and large amplitude oscillation states are localized in the 
right and left-hand wells, respectively, and are separated 
by a "potential barrier" with purely imaginary momen- 
tum. Here we do not observe any interference effects, as 
opposed to the case of the parametrically driven oscilla- 
tor [8j . We apply the WKB theory in this "potential" in 
order to determine the tunnel splitting in the limit of a 
low transmission through the forbidden region. The clas- 
sical turning points Xi are given by ps,L{xi,E) — 0, see 
Fig.[l](a). The bound state energies at zero transmission 
are given by the Sommerfeld energy quantization rules 

Si2{E) = 7rn+7r/2, Si-z'iE) = 7rm+7r/2, n,rn G (5) 

where Sij{E) — J^' sign{x — X)\p{x, E)\dx/h and the 
negative sign on the left hand side of Xis due to mir- 
roring. Whenever a pair of energies from either well is 
degenerate, resonant tunneling through the barrier can 
occur. This induces coupling between the two wells and 
lifts the degeneracy. The level crossings become avoided 
crossings at finite transmission and the full WKB condi- 
tion reads 

cotSi2{E)cotSi,3'{E) = exp{-2S3'iiE))/4. (6) 

We expand the quasi-energy E and the actions Sij in 
series of ^ = l/4exp (— 2S'3/i) around the level crossings 
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FIG. 2: (a): Quantized energies: eig envalues (EV) of H^*' 
versus WKB. H^"'> was represented in the number state ba- 
sis considering a number of levels 2N. (b): Tunneling- 
induced energy splittings at level crossings. Frequency sweep 
at mfl/h = 2, 7 = mQ,^/24, kloc/O.^ = 0.1 and Fo = 0.5Fc(//). 



with quasi-energy Eq where eqs. ([5| are simultaneously 
satisfied. The first energy correction Ei£^ is obtained 
straightforwardly from dESi2\EodES4'3'\Eo{Ei£,)'^ = 
and the tunneling rate is obtained directly from the en- 
ergy splitting at the avoided level crossings 



Tt = 



exp(-53'i) 



(7) 



This can be evaluated in closed form involving elliptic 
integrals for Sij and we obtain the exact expressions 
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(^)m wWi and e e {m,s,M}. 
Thus, for at (m) and S'3/1 at (s) we reproduce the 
harmonic oscillator result. The saddle point "frequency" 
Os is imaginary as expected. 

We simphfy Eq. (H) by locally approximating h'^^ 
close to the extremal points by harmonic oscillators, 
i.e. assuming that Sij are linear functions of E. This ap- 
proximation holds for all Sij simultaneously when E is far 
enough from both extremal points -Es.m, as it is the case 
for the ground state E^ + of the small amplitude 

well. In this approximation Syi{E) k, ti{Es — E) / {h\Q.s\) 
and thus we find a compact approximation 

Ft « f^„,/7r2exp(-^(S, -E^- hnj2)/{h\n,\)). (8) 

Our calculations rely on a series of assumptions. To 
test them, we compare the results to a full numerical 
diagonalization of h'^^ taking a basis of the first 2N Fock 
states. At = 0, the number of levels that cover the 
bistability region is = hil{2mil)^ / . As shown 
for a representative set of data in Fig. [2] we find good 



agreement between these numerically exact results and 
the predictions of Eqs. (5[7) and also ([s]). 

Quantum tunneling is significant only close to level 
crossings. It always competes with the activation over 
the barrier, which occurs at all energies and is based on 
classical fluctuations due to coupling to a heat bath. A 
rather detailed treatment of a similar process has been 
given in Refs. [10]. We now estimate these effects and 
compare them to the quantum tunneling rate. When 
modeling activation, it is crucial to consider that we are 
working in a frame rotating relative to the heat bath, 
which is fixed in the laboratory. 

We start from Eq. ([2]) . As we will adopt the mean first 
passage time approach |16j . it is sufficient to approxi- 
mate the system Hamiltonian close to its minimum in 
phase space by H^f^ ~ p'^ / {2mcs) + V (x) where the effec- 
tive mass is determined by the curvature of the Hamilton 
function m 
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(a;,p)|m and the effective poten- 
tial is V{x) = Hq"' (x^pni). In this approximation we 
obtain a quantum Langevin equation 

f°° 2J(uj) /"* 
mcf[X+dxV{x)—x du! K^ricft / k{t—s)x{s)ds — £,{t) , 



where 

k{t) = 



2J{uj) cos((tJ — i')t) 



(cj - v)nmcs 



duj. 
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k{t) is peaked on a short time scale lo^ ^. Its magnitude 
is characterized through the effective friction constant 
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k{t)dt = 2k (5 
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The factor of two difference between Koff and the damp- 
ing constant of the undriven harmonic system accounts 
for the fact that in the rotating frame there are bath 
modes above and below = (see Eq. ([2|) whereas 
for the undriven case the frequencies are strictly posi- 
tive. Thus oscillators with frequency uj have the spectral 
density J{lo + v) and modes with negative frequencies 
have significant contribution to noise even at low tem- 
peratures. We use a detailed balance condition to de- 
termine the effective temperature of the bath as seen by 
a detector in the rotating frame, e.g. a two level system 
with level separation Ml-^ 



p{n^, T)/p{-n^, T) = exp {hn^i3,s) ■ 



(9) 



Here P{uj, T) = J{uj + + n{ui + v, T)) is the proba- 
bility for a quantum Tiio to be emitted to the bath in the 
rotating frame. The increase of the quasi-energy in the 
rotating frame can actually lower the energy in the lab- 
oratory frame, i.e. what a detector in the rotating frame 
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regards as absorption can actually be emission in the lab 
frame. In the case of constant acceleration in relativistic 
context this behavior is known as the Unruh effect fTfl. 
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FIG. 3: (a) The ratio of tunneling and activation rates from 
the small well at the avoided level crossings, (b) Correspond- 
ing tunneling rates compared to KcEF (where Ft /Fa > 1). Driv- 
ing frequency sweep at Fo = O-YFciv); Values of Q (in GHz): 
!(•), 2(n), 3(»), 4(a), 5(t), 6(>), 7(+), 8(x), at parameters 
specified in text. 

The barrier crossing problem for systems described by 
a quantum Langevin equation is well studied in the con- 
text of chemical reactions. For low damping, ^ ilai 
mean-first-passage time theory predics the activation rate 
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where S{E) — § p{x,E)dx. In the traditional low tem- 
perature limit KcsS{Es) <C ksTcs <C £^s ~ -^'m the acti- 
vation rate becomes 

Ta = Keff/3efff^mexp (-(£;, - E^)p,s) S{E,)/{2tt). (11) 

In our case, the noise temperature ksTcn can be larger 
than the barrier height Eg — In this limit we obtain 



from Eq. (lOl 



Fa = Keff (F(/3eff (^« " ^^m)))" 



(12) 



where F{x) = J (ix(exp(a;) — l)/x = Ei(a;) — log(a::). 

Summarizing, in the rotating frame, as a consequence 
of driving, the bath appears with a quality factor fim/Koff 
reduced by approximatively a factor of two and an en- 
hanced effective temperature Tog. Moreover, the bath 
shifts the detuning 5. We show that experimental obser- 
vation of MDQT could still be possible. At the the level 
anticrossings we calculate the WKB tunneling rate from 



the ground state and the activation rate from Eq. (12 1, 
see Fig. [3]^ a) where we have considered a Josephson junc- 
tion with K — 10^^f2, the temperature T = 10 mK, shunt 
capacitance C = 2-10^^^ F and 7 = TOri^/24. The values 
of 5 where these anticrossings occur are found by min- 
imizing I cot(S'4'3'(i?in))| and are in agreement with the 
weak driving result [H] i5 = 37n/(2m^ri^), n G N. We 
observe that the quantum tunneling rate can be one or- 
der of magnitude larger than the activation rate in the 



limit of relatively small detuning 5 and low damping. By 
increasing the value of a = mfi /fi, we observe a reduction 
of the ratio Ft /Fa as expected, since a measures the num- 
ber of quantized levels in the system and thus the "classi- 
cality" of its behavior. In Fig. 3]we have a S (2, 20), while 
in the experiment of Ref. [13] a was larger than 100, at 
higher temperature and smaller quality factor, such that 
MDQT was probably masked by thermal activation. We 
expect that at the values of Fig. [3] the experiment we 
propose should produce direct evidence for MDQT. 

In conclusion we have investigated macroscopic dy- 
namical tunneling by mapping it onto tunneling between 
two potential surfaces and found a compact analytic ex- 
pression for the tunneling rate. We compared this process 
with the activation over the barrier using the mean first 
passage time approach. The values obtained suggest that 
dynamical tunneling can be singled out from the back- 
ground of activation processes. We have proposed an ex- 
periment realizable within existing technology to demon- 
strate dynamical tunneling by monitoring the switching 
rate between the two dynamical states while tuning a 
parameter of the external driving. 
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